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Abstract 

We model the disassembly of an excited nuclear system formed as a result 
of a heavy ion collision. We find that, as the beam energy in central colli- 
sions in varied, the dissociating system crosses a liquid-gas coexistence curve, 
resulting in a first-order phase transition. Accessible experimental signatures 
are identified: a peak in specific heat, a power-law yield for composites, and 

a maximum in the second moment of the yield distribution. 
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Nuclear matter is a fictitious arbitrarily large N = Z system in which the Coulomb 
interaction is switched off. Mean-field theory for nuclear matter has been done many times, 
and it is well-known that such a system shows a Van-der-Waals-type liquid-gas phase tran- 
sition. It was suggested in the early eighties that in heavy ion collisions at intermediate 
energies one might be able to probe this liquid-gas phase transition region. In heavy ion 
collisions, matter would be heated as well as compressed. This compressed blob would then 
expand passing through the liquid-gas coexistence phase. One might be able to extract in- 
formation about this region from selected experimental data. Earlier, the Purdue group had 
conjectured that the break up of large nuclei by energetic protons would show signatures of 
critical phenomena 0. 

There are several complications which make the study of phase transitions in nuclei 
difficult. The systems are small, thus singularities get replaced by broad peaks. The collision 
is over quickly. The existence of thermal equilibrium has sometimes been questioned and 
replaced by quite complicated transport equation approaches. Often one has been content 
to do calculations where the objective is to fit the experimental data. Once this is achieved 
the question of a possible liquid-gas phase transition is not addressed. There are models 
where such questions are irrelevant, or at least very indirect, such as various models based on 
sequential decays. The literature on calculations for fragment yields in intermediate energy 
heavy ion collisions is huge. We will not attempt to mention all approaches. 

In this paper we will focus on the liquid-gas phase transition using a lattice gas model 
0. Previously we have used the model to fit data on central collisions 0, on peripheral 
collisions |Q, and for central Au on Au collisions ||^. The model gives a fair description of 
data in those instances — does it say anything definite about the experimental signature of 
liquid-gas phase transitions? 

In the lattice gas model we place n nucleons in cubes, where n is the number of nucleons 
in the disassembling system, and N/n = po/Pf, where po is the normal nuclear density, and 
the disassembly is to be calculated at pf, the "freeze-out" density beyond which nucleons are 
too far apart to interact. An attractive nearest-neighbor interaction is assumed. We place 
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the n nucleons in N cubes by Monte-Carlo sampling using the Metropolis algorithm. Once 
the nucleons have been placed we also ascribe to each of them a momentum. The momenta 
are generated by a Monte-Carlo sampling of a Maxwell-Boltzmann distribution. Various 
observables can be calculated in a straightforward fashion. Experiments usually measure 
the distributions of cluster sizes, i.e., the yield Y{Z), which is a function of the number 
of protons in the composite Z. In our model, two neighboring nucleons are considered to 
be part of the same cluster if their relative kinetic energy is insufficient to overcome the 
attractive bond: -|- e < 0. This can be immediately turned into a temperature and e 

dependent bonding probability much like Coniglio-Klein's [Q. The prescription allows 
us to calculate the distribution of clusters. 

The equation of state in the lattice gas model has been well studied in condensed matter 
physics. The grand canonical ensemble for the lattice gas corresponds to a three-dimensional 
Ising model in the presence of a magnetic field. It is thereby possible to translate many well 
known results to the present situation. Fig. 1 depicts the phase diagram of the lattice- 
gas model in the thermodynamic limit, which is adapted from Ref. In drawing the 
coexistence curve DCE, the series expansion given in Ref. |0] was used. The point "C" is 
the thermal critical point which occurs at T 1.1275|e| and density p//po = 1/2- The 
coexistence curve is DCE, and CB is the line along which percolation sets in. The line CB 
is only slightly different in our and Coniglio-Klein's prescription. 

In an experimental situation involving the collision of heavy ions, one gates on central 
collisions and varies the beam energy. At low beam energy, the freeze-out density will 
be below the co-existence curve. As the beam energy increases, this point will cross the 
coexistence curve (indicated by an arrow in Fig. 1). Note that there is no reason for the 
system, under heating, to "tune itself" to the second-order point at density p = 1/2, since 
there is no symmetry between the clustered and unclustered phases involved in the heavy-ion 
collision. Indeed, in statistical models of disassembly of which we know PJTO[], the freeze-out 
density is less than O.Spo- In our model, we find that the data are best fitted by a pf between 
0.3po and 0.4po- Therefore the phase transition is first order. 
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As one crosses the coexistence curve, experiments would see various signatures of the 
first-order transition. In the thermodynamic hmit these are well defined. For example, 
since first-derivatives of the free energy are discontinuous in the thermodynamic limit, there 
are delta-function peaks in second derivatives |Tl||T^, e.g., in the heat capacity at constant 
pressure Cp, and in the isothermal compressibility k^. There is also a clear signature of the 
transition in the density- density fiuctuation correlation function r(r) (whose integral gives 
the compressibility by a thermodynamic sum rule |jl3| / d-r r{r) oc kt)- This correlation 



function decays exponentially above and below the transition point. At the transition point 
itself, however, r(r) is fiat |T^JT^ , as is evident from, for example, the sum rule. 

This singular behavior is smeared out in finite-size systems [|14[. Instead of delta- function 



singularities, kt and Cp have broad Gaussian peaks at the first order transition, whose 
height is proportional to the system's volume ~ L^. Indeed, these features are somewhat 
analogous to those at a continuous second-order transition: In the thermodynamic limit at 
a continuous transition, Cp and kt have power-law singularities, but in a finite-size system, 
these singularities are replaced by bumps of height ~ L^^'^, where 7 and u are critical 
exponents. This similarity implies why it is often difficult to distinguish the two types of 
transitions in a finite-size system [|r^,Q. 

While Cp and kt are both infinite at a first-order liquid-gas transition, the heat capacity 
at constant volume C^ has only a step discontinuity in the thermodynamic limit, giving a 
bump in a finite-size system. The extent of the discontinuity is given in terms of thermo- 



dynamic relations in Ref. [|T^. This is in contrast to the more pronounced behavior at a 
second-order transition where C^ ~ L°'^'^, and a is a critical exponent. 

In Fig. 2 we show results of our numerical calculation. We fix at 7^ and vary n to 
obtain a variable freeze-out density. Half of the nucleons are labeled as protons. Although 
this system is small, there is a well defined bump in C^ signifying the transition. An arrow 
marks the nearby point at which the transition occurs in the thermodynamic limit. 

In nuclear experiments locating the peak in Cy as the beam energy increases is difficult, 
although in recent years tremendous progress has been made in the measurement of temper- 



ature |T^. However, other characteristics of the transition appear which are more readily 
measurable. In particular, the yield Y{Z) is readily measurable, which gives the distribution 
of clusters of charge Z, n{Z). The probability of clusters of a given size is related (through 



its second moment), to the density correlation function p|, p!7|JT8[] : Roughly, if the clusters 
are not fractals, Z'^n{Z) ~ r(r), in the disordered phase, where r is the diameter of a cluster 
of charge Z. At the transition, a droplet of thermodynamic size spans the system. Hence 
the yield includes an infinite cluster, or in percolation language, a cluster which spans the 
system. The remainder of the distribution describes fluctuations in that phase, giving the 
density correlation function, whose integral is the compressibility. The delta-function peak 
in the compressibility implies that the correlation function, and hence the yield, is broad. 
Following standard practice [|19|, the yield is fit to a power law form Y{Z) oc 1/2"^, giving 



an effective exponent r, even when the distribution has deviated from a power law. In fact, 
at a continuous phase transition, t = 2 + f3/{(3 + 'j), where /? is a critical exponent. Here, 
however, where the transition is first-order, the correlation function decays exponentially in 
either bulk phase (i.e., the effective r ^ oo above or below the transition), while, since the 
correlation function is fiat at the transition itself, and is related to the second moment of the 
cluster distribution, r = 2 there, for an infinite system. Of course, in a finite-size system, 
there will be an effective r, which is neither two nor infinity, so that the transition will look 
somewhat analogous to a continuous transition. 

It is also useful to consider the second moment of the cluster distribution function, 
5*2 = J2'yiA'^n{A)/n, where n{A) is the number of nucleons with mass number A, and n is 
the total number of nucleons. The primed sum excludes the largest cluster, and so we expect 
it to be proportional to the compressibility kt- The usefulness of the second moment 
was emphasized by Campi pO |. 



Fig. 2 shows that these quantities do indeed provide clear experimental signatures of the 
first-order transition in a finite-size system. The crossing of the coexistence line is evident 
in the minimum value of the effective r, which is close to 2, and the prominent maximum 
in 5*2 and Ct,. All of these occur at approximately the same point, close to the transition in 



the thermodynamic hmit, indicated by an arrow. 

For our analysis, it is important that the freeze-out density is on the low density side of 
the coexistence curve. Beyond the critical point, "C", the bump in C^, continues on the line 
CE whereas the minimum of r will follow line CB, which corresponds to a change in short- 
range order. We should also note that very large Coulomb forces will alter this picture. We 
now describe briefly those effects, following Ref. In our lattice gas model one generates 
an event in which the nucleons have been placed with appropriate momenta. Using our 
prescription we can immediately obtain the cluster distribution. Starting from this config- 
uration we can propagate the system by molecular dynamics with a suitably chosen short 
range nucleon-nucleon interaction (to correspond to the nearest neighbor interaction) with 
or without the Coulomb interaction. By propagating molecular dynamics for a considerable 
time one can unambiguously obtain the cluster distribution, and determine the effect of the 
Coulomb interaction. Ref. concluded that without the Coulomb interaction there is little 
difference between lattice gas results and molecular dynamics results, but that Coulomb in- 
teractions can be important for large mass numbers. For mass number A = 85, the Coulomb 
interaction is a small perturbation, while a.t A = 394 it is so strong as to totally alter the 
picture. In the latter case there is no minimum in r; at a temperature of 1 MeV, the value 
of T is slightly higher than 1 and it continues to rise monotonically with temperature. We 
have now studied this in much greater detail and find that a minimum in the value of r 
continues to be obtained for mass number as large as A = 200. Details of this as well as 
further applications of the lattice gas model will be published in a longer paper. 

To conclude, we have modeled the disassembly of nuclear matter following a heavy ion 
collision. We find that the transition is first order, with the standard signatures of such a 
phase transition. 

We thank Charles Gale, Martin Zuckermann, and Hong Guo for useful discussions. This 
work is supported in part by the Natural Sciences and Engineering Research Council of 
Canada and by le Fonds pour la Formation de Chercheurs et I Aide a la Recherche du 
Quebec. 
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FIGURES 

FIG. 1. Phase diagram of three dimensional lattice gas model. The line DCE is the coexistence 
curve. CB is the percolation line. The arrow demonstrates the crossing of the coexistence line as 
the beam energy for central collisions increases. In the Ising model DCE is the line of spontaneous 
magnetization. With the usual convention the point on the abscissa corresponds to magnetization 
-M and the point 1 to M. 

FIG. 2. Curves for r, Cy, and 52 at Pf/po = 0.2, 0.3, and 0.4. At each point 1000 events 
were taken. is in units of the kinetic energy contributes 1.5 to it at all temperatures; 
Tc ^ 1.1275|e|. 
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